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We investigate electron interferometry of edge states in Topological Insulators. We show that, 
when inter-boundary coupling is induced at two quantum point contacts of a four terminal setup, 
both Fabry-Perot-like and Aharonov-Bohm-like loop processes arise. These underlying interference 
effects lead to a full electrically controllable system, where the magnitude of charge and spin linear 
conductances can be tuned by gate voltages, without applying magnetic fields. In particular we 
find that, under appropriate conditions, inter-boundary coupling can lead to negative values of the 
conductance. Furthermore, the setup also allows to selectively generate pure charge or pure spin 
currents, by choosing the voltage bias configuration. 
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I. INTRODUCTION 

Interferometry is a hallmark of wave physics. While 
optical interferometers have been known since past cen- 
turies, more recently remarkable efforts have been made 
to realize electronic interferometers, which exploit the 
wavelike nature of electrons, a signature of quantum 
world. Due to the decoherence effects that electrons 
experience in ordinary macroscopic metals, such type 
of wavelike phenomenon can be observed only in meso- 
scopic systems. Fabry-Perot interference pattern caused 
by electron waves reflected at the contacts between a 
carbon nanostructure and two metallic electrodes have 
been observed in various experiments on nanotubeiP^ 
and graphene^. Similarly, clear evidence of Aharonov- 
Bohm effect has been demonstrated in semiconductor 
heterostructures and graphene rings threaded by a mag- 
netic fluxPEI and various electron interferometers have 
been realized in Quantum Hall Effect (QHE) systemsPE] 

Electron interferometry is not just a conceptually im- 
portant effect, it has significant practical consequences, 
for it can be exploited to realize quantum transistors 
with high current-carrying capability^. Indeed currents 
can be switched 'on' and 'off' by varying the electron 
interference conditions from constructive to destructive 
through a gate voltage or a magnetic field. 

At present, most electronic interferometers are based 
on the electron charge. However, electron is also charac- 
terized by its spin, and spintronics - the field investigat- 
ing transport and manipulation of information with such 
degree of freedom - is experiencing an extremely rapid 
growth. Indeed spin is much more robust to environment 
decoherence effects with respect to electron charge^, and 
spin coherence lengths may reach and exceed lOO/miP^l 
For these reasons, interferometry involving spin has been 
proposed to realize transistors and filters exploiting spin- 
orbit induced precession in semiconductord^Ell, ferro- 
magnetic materials^HUl anc l magnetic fields^. In spite 
of these advances, interferometric control of spin cur- 
rents remains a difficult task requiring ad hoc optical 



technique s] 25 * 26 !, and has not reached the state of the art 
level. Thus, as far as interferometry is concerned, spin- 
tronics is not as competitive as ordinary charge-based 
electronics yet. 

A major boost to spintronics is expected to come 
from the recent discovery of Topological Insulators. A 
Topological Insulator (TI) is a bulk gapped material ex- 
hibiting conducting gapless channels at the boundaries.^ 
In these edge states the generation of spin currents is 
greatly facilitated from the close connection between mo- 
tion direction and spin orientation (helicity). In two- 
dimensional realizations of a TI, for instance, only spin-f 
electrons propagate rightwards and only spin- 4- electrons 
leftwards, along a given boundary. Remarkably, TI edge 
states behave as perfectly conducting one-dimensional 
ballistic channels, since impurity backscattering is pre- 
vented from time reversal sym metry. These peculiar 
properties, theoretically predicte d 28 * 29 * 31 ! and experimen- 
tally observed in HgTe/CdTe quantum wells' and in 
various ot her mate rials^, make TIs ideal candidates for 
spintronics ) 19 l 27 l 34 l Ml 

For these reasons the investigation of electron inter- 
ference effects in these systems appears to be a particu- 
larly timely issue. Quite recently, for instance, evidence 
of Aharonov-Bohm interference pattern in magnetoresis- 
tance of TIs has be en exp erimentally observed^! and the- 
oretically discusse d 35 * 38 ' . Notably, comparative analysis 
have pointed out that, with respect to the cases of edge 
states in QHE bara^l and SU(2) symmetric systems-^, 
interferometry of TI edge states exhibits intrinsically dif- 
ferent behavior, and represents a challenging open prob- 
lem with its own peculiarities. 

So far, most studies concerning TI in terfer ometry have 
involved magnetic flux or ferromagnets I25133ES1 However, 
TI edge states appear even in the absence of magnetic 
field, since the mechanism underlying their properties is 
the spin-orbit coupling and the related inversion of the 
electronic band order. Indeed this represents one of the 
crucial advantages of TIs with respect to QHE edge states 
in view of device miniaturization, due to the difficulty in 
realizing scalable circuits operating under high magnetic 
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fields or space varying magnetizations with nanoscale res- 
olution. It is thus desirable to explore TI interferometry 
also in absence of magnetic fields or ferromagnets. 

In this article we address interferometry based on 
purely electrical effects, i.e. in the presence of time rever- 
sal symmetry. We shall focus on a two-dimensional real- 
ization of a TI, and consider the interferometer sketched 
in FigJT] where Kramers pairs of TI edge states flow 
at the Top and Bottom boundaries of a quantum well. 
The electrochemical potential of each injected edge state 
is controlled by a metallic electrode, characterized by a 
voltage bias Vi (i = 1, ...4). Inter-boundary scattering 
between edge states can occur at two Quantum Point 
Contacts (QPCs), giving rise to the possibility of loop 
trajectories, which determine the currents through elec- 
tron wave interference. Here the interference conditions 
are set by two (Top and Bottom) gate voltages V a> T and 
Vg t B, which modify the electron phase in the loop pro- 
cesses by shifting the electron momenta in the regions 
between the two QPCs. The setup of Fig[l]can thus lead 
to a full electrically controllable system, where charge 
and spin conductances of each electrode can be tuned by 
the gates. We shall determine and discuss the behavior of 
the conductances as a function of V 3i t and Vg,B, for vari- 
ous configurations of the biases Vi. In particular we shall 
show that, under appropriate conditions, inter-boundary 
coupling can even lead to negative conductance values. 
Furthermore, the setup also allows to selectively generate 
and tune pure charge or pure spin currents. 

The article is organized as follows. In Section [TT] we 
present the model for the setup, in Section [TTT] we discuss 
the interference phenomena of the system and point out 
the relations with other electron interferometers. Then, 
in Section |TV| we show the results concerning charge and 
spin currents, mainly focussing on two specific configu- 
rations of the four terminal setup. Finally, we discuss 
the results in Section Ivl and draw our conclusions in Sec- 
tion rvn 



II. THE MODEL 

The edge states at the Top and Bottom boundaries of 
the device are characterized, at low energies, by a linear 
spectrum described by the following HamiltoniarP3 

H = -ihV F ^ J dx : ^Rai x ) d ^ Rai x ) ■ 

-:^(x)d x *^(x) :1 (I) 



<r=t,r 



where x denotes the longitudinal coordinate, and 
\& La the Right and Left mover electron field operators, 
and er =y,i the spin component. Here we adopt the no- 
tation that spin-f Right (Left) movers and spin-J, Left 
(Right) movers flow along the Top (Bottom) boundary, 
as depicted in Fig [I] Without loss of generality we shall 
assume that the equilibrium Fermi level Ep of the de- 
vice in the absence of any bias is located at the Dirac 




FIG. 1. (Color on line) Schematic description of the pro- 
posed four terminal setup, where edge states flow at the Top 
and Bottom boundaries of a TI quantum well. Blue (red) 
lines denote spin-t (\) edge channels, and are slightly sep- 
arated for illustrative reasons. The dashed box denotes the 
scattering region, where inter-boundary tunneling occurs at 
two QPCs, located around positions x a and Xh and separated 
by a distance L = x b — x a . Two gate voltages V 3 ,t and V 9 ,b 
shift the edge state momenta in the Top and Bottom regions 
between the QPCs, modifying the electron phase in the loop 
processes induced by tunneling. No magnetic flux is present. 



point of the spectrum. The symbol : : in Eq.Q denotes 
the normal ordering with respect to the equilibrium state 
where all levels below Ep axe occupied. Variations from 
this energy level can be induced by the gate voltages, as 
discussed below. 

The constrictions of the QPCs induce inter-boundary 
scattering. It can be shown with quite general 



argument; 
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that time-reversal symmetry only allows 



two types of tunneling terms, namely a spin-preserving 
tunneling 

«Ln= E Jdx(T p (x)^ Ra (x)* LlT (x) + 

+r;(x)*l-(x)* m (x)) (2) 



and a spin-flipping tunneling 



H f tun = £ a J dx ( r /( a ) + 



+r* f (x)n>l ti (x)y ctt (x)) (3) 



a=B/L=± 



In fact, the pinching of the two edges states at the QPCs 
causes a local modification of the spin-orbit coupling with 
respect to the bulk case , so that both terms are ex- 
pected to contributePEH In Eqs.Q and Q T pJ (x) de- 
note space-dependent tunneling amplitude profiles. The 
two QPCs will thus be described through a profile peaked 
around two centers x a and Xf, (see Fig{T]), and rapidly 
decaying beyond a longitudinal lengthscale £. We shall 
assume that the constriction are short with respect to the 
Fermi wavelength Xp, and that the distance L = Xf, — x a 
between the two QPCs is large compared to £, i.e. £ < 
\p < L. Under these conditions, it is sufficient to assume 
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that tunneling is point-like. For the sake of simplicity we 
shall also consider that the two QPCs are characterized 
by equal tunneling amplitudes, so that the profiles can 
be taken as 



(4) 



where j p and 7/ are real dimcnsionlcss spin-preserving 
and spin-flipping tunneling amplitudes, respectively. 

Finally, the coupling to the two gate voltages V g _r 
and V g B can be described by the term 



dx [ eV g .T ( p m {x) + p Li (x)) + 

+ eV 9,B (Pri(x) +p Lt {x))] 



(5) 



where 



P a a=^L{x)^ ar7 (x) 



(G) 



denotes the electron density, with a — R/L and a =t,4- 
As anticipated above, V 3 ,t and V 9t B shift the electronic 
spectrum ([I]), and their difference breaks the degeneracy 
between Top and Bottom boundaries. 

The equations of motion for the four fields f acr are 
easily obtained from the total Hamiltonian of the device 



% — V.Q + %tun + ^tun + Hg (?) 



where the four terms are given in Eqs.Q, @, ^ and (§, 
respectively. For ^ R ^, for instance, one obtains 

d t ^ m = v F (- d x ^ R -f - ikg t Rf9(x a < x < x b ) ^ m (x)- 
-2i (1p *Lf(as) + If *m,(x))S(x - xi) (8) 

l=a,b J 

where 9 is the Heaviside function. Similar equations are 
obtained for the other fields. The four equations, coupled 
by the inter-boundary tunneling terms, can be solved by 
superposing solutions corresponding to fixed energy val- 
ues. At a given energy E, measured with respect to the 
equilibrium Fermi level Ep , the solution can be obtained 
by the Ansatz 
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where k E = E/hv F , k 3>R ^ = fc 9 ,ij. = eV g , T /ftv F , k g ^ = 
k g ,Lt = eVg^B/fcvp, and cp a a denotes the electron op- 
erator inside the scattering region. The incoming and 



outgoing electrons, respectively described by the opera- 
tors a and b in (|9j), are connected via the entries Sij of 
the Scattering Matrix JpU 
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In the columns of Eq.(10l a and b operators appear ac- 
cording to the clockwise order of the four leads in the 
setup Fig[T] A lengthy but straightforward calculation 
allows to determine the Scattering Matrix entries Sij as 
a function of the tunneling amplitudes 7 P ,7/ and the 
gate voltages V 3 ,t an d V 3t B- Helicity and time-reversal 
symmetry lead to vanishing diagonal entries, 



Su = 



(11) 



while the other entries can be expressed in a compact 
form by utilizing the transmission coefficients of each sin- 
gle QPC in terms of the tunneling amplitudes j p and 7/, 



rpa 

1 2l 



1 3l 



rpa 

1 34 



rpa 

1 42 



rpa rpa 

1 4l — 1 32 





(1 












(1 






(1 


'P 




(1 




1}? 



(12) 



Here T, a 4 (T/i) denotes the transmission coefficient from 
lead j to lead i of the Left (Right) QPC alone. Time 
reversal symmetry ensures 



rpl rpl 



I — a,b 



(13) 



Furthermore, we have assumed that the QPCs have iden- 
tical parameters 



rpQ 

ij 



rpb 

ij 



(14) 



Then, the S-matrix entries read 



S31 — 



1 + Vt^I e 2 ^--^^ 



S32 — 

S34 — 
where 



1 + V T 43 T 2i e MkE ~ 



-1 e 



-2ikEX} ) 



ppb~ 1 /rpa „2i{k E 
1 34 + V J 34 e 



(15) 



-)L 



1 + V T & T 2X eMkE ~ 



g,T/i 



hvi 



(16) 
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Similar expressions are found for the other entries, to- 
gether with the following relations 



\Si 3 (E)\ = \S 31 (E)\ = \S 24 (E)\ = \S 42 (E)\ 



\Su(E)\ = \S 41 (E)\ = \S 23 (E)\ = \S 32 (E)\ 



(17) 



The S matrix entries allow to define the system trans- 
mission coefficient from lead j to lead i (computed at the 
Fermi energy, i.e. at E = 0) as 



\Si 



(18) 



It is well knowrl^l that Tij appear in the expressions for 
the current, as will be explicitly discussed below. 



III. INTERFERENCE PHENOMENA 
IN THE SYSTEM 

Before presenting the results for the currents (see Sec- 
tion 



IV I, in this section we describe the interference ef- 



fects characterizing the setup. This gives us the opportu- 
nity to point out similarities and differences with respect 
to other electron interferometers. 



A. Fabry-Perot interference 

To simplify the discussion let us first assume that 
the spin-flipping tunneling amplitude is vanishing 
(7/ = 0). In this case quantum interference is due 
to spin-preserving processes 0y p 7^ 0). An illustrative 
example is described in Figllja): a spin-j" electron is 
injected from terminal 2 and, when reaching the right 
QPC, partly tunnels to the Bottom boundary as a 
left-mover; it then tunnels back to the Top boundary at 
the left QPC and eventually reaches terminal 3. The 
phase accumulated in the loop causes the interference 
with the electron wave that straightforwardly travels 
from terminal 2 to 3. Similar processes occur for 
spin-J, electrons injected from terminal 1 to 4. Notice 
that these processes involve a change in the motion 
direction, i.e. an energy dependent momentum transfer 



ikp-x 



-ik E x^ with ks _ E / hvp _ 



This effect is reminiscent of the Fabry-Perot (FP) i n- 
terference pattern occurring in carbon nanotubes^HH r 
in single channel quantum wired^, and involving 2k p mo- 
mentum transfer. Indeed the inter-boundary terms ^ 
causing the loop of Fig{2^a) correspond to backscatter- 
ing terms originating from the contact resistance at the 
nanotube(wire)/electrode interfaces, while the distance 
L between the two QPCs plays the role of the length of 
the nanotube (wire). 

The electron phase accumulated in the loop depends 
on the gate voltages, which modify the momenta in the 



Top and Bottom region between the QPCs. In view of 
the analogy discussed above, we shall denote such phase 



>fp = 



e(V g ,T + V 9 ,b)L 
hvp 



(19) 



as the Fabry-Perot (FP) phase of the device. Notice 
that 4>fp depends on the sum V 3t B + Vg,T of the gate 
voltages and on the QPC distance L. 



B. Aharonov-Bohm interference 

Let us now consider the effect of the spin- flipping terms 
(7/ 7^ 0), and assume that the spin-preserving processes 
are absent (j p = 0). Interference processes for this sit- 
uation are schematically depicted in Figj2^b) : A right- 
moving spin- 1 electron injected from terminal 1 exhibits 
two possible paths to reach terminal 3, corresponding to 
tunneling to the Top boundary occurring at the QPC on 
the left and on the right. The phase difference 



'AB 



e(V g ,T - V 9 ,b)L 
Hvf 



(20) 



between these two possible paths connecting terminals 1 
and 3 causes an electron interference phenomenon similar 
to the electrostatic Aharonov-Bohm (AB) eff ect o ccur- 
ring in semiconductor and metallic rings rings pOSl Sim- 
ilar AB-like loop trajectories connect terminals 2 and 3. 
Notice that, differently from the spin-preserving FP loop 
processes, the AB interference due to spin-flipping pro- 
cesses is independent of the energy E and involves the 
difference V s ,T—V g .B between the two gate voltages. This 
is due to the fact that, while in FP processes the elec- 
tron travels along the two arms of the loop in opposite 
directions [see Figj2^a)], for spin-flipping processes the 
electron always preserves its motion direction in order for 
the spin to flip [see Figj2jb)]. This constraint originates 
from time reversal symmetry, which leads spin-preserving 
inter-boundary tunneling to occur backwards [R <-> L, 
see Eq.Q], and spin-flipping inter-boundary tunneling 
to occur forwards [a — > a, see Eq.(|3])]. This connection 
between spin orientation and motion direction is a hall- 
mark of topological insulator edge states dynamics. 



C. General case 

The two processes of Fig[2] described above refer to 
the extreme cases where either "f p or 7^ vanishes. In the 
general case (jp 7^ and 7/ 7^ 0) higher order processes 
give rise to multiple loops where the Fabry-Perot and the 
Aharonov-Bohm interference effects interplay with each 
other. These processes directly enter the Scattering Ma- 
trix entries ( 15 ) and thus the transmission coefficients 



(18). In particular, to lowest order in the tunneling am- 
plitudes 7 P and 7/, the 'horizontal' transmission T 32 is 
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that the gate coupling term ^ can be rewritten in the 
following way 




FIG. 2. The two types of basic interference processes allowed 
by time reversal symmetry of the system, (a) An example of 
loop process induced by the spin-preserving tunneling terms 
Eq. pj, reminiscent of a Fabry-Perot interferometer. The 
electron flows rightwards in the Top Boundary and leftwards 
in the Bottom boundary, so that the loop interference phase 
is <j) FP = e(V g ,T + V g , B )L/hv F . 

(b) An example of loop processes induced by the spin- 
flipping tunneling terms Eq. pj, reminiscent of electrostatic 
Aharonov-Bohm interferometer. The electron maintains the 
motion direction both in the Top Boundary and in the Bot- 
tom boundary, so that the loop interference phase is 4>ab = 
e (Yg,T — V g ,B)L/hvF- A similar loop connects terminal 2 to 3. 
In general multi-loop processes, where the two processes in- 
terplay, are also possible. 



determined by both Fabry-Perot and Aharonov-Bohm in- 
terference processes, whereas the 'crossing' transmission 
T31 is only affected by Aharonov-Bohm loops. Indeed, 
while terminal 2 can be connected to terminal 3 through 
either one FP loop or one AB loop (see Figj2]), terminals 
1 and 3 can only by connected through one AB loop, 
as can be formally checked by perturbative expansion of 
Eqs.(fl5j) and (pi. 



D. Differences from other electron interferometers. 

In the previous subsections we have highlighted the 
analogies between the basic interference phenomena of 
the setup Figjl] and other electron interferometers. Here 
we wish to point out the differences. 

The first difference emerges at the level of the cou- 
pling to the gate. To illustrate this point, one can realize 



dx[(V gtT + V g>B )p c (x)+ (21) 
+ v F \V g , T -V g , B ) I s (x)] 



where 



p c (x) = e (pat(x) + PLt(x) + Pr±{x) + PlU x )) ( 22 ) 

is the charge density operator, and 

I s (x) = ev F (pRf(x) - pl T (x) - Pm(x) + Pli(x)) (23) 

is the spin current operatoi^, with the electron densities 
Pan given by Eq.([6]). Equation (21) shows that the 



sum of the gate voltages, i.e. the FP phase (19 1, 
couples to the charge density, while their difference, 
i.e. the AB phase (20), couples to the spin current. 



The expression (21) points out a difference with respect 
to Fabry-Perot electron interferometers realized with 
carbon nanotubes or quantum wires. In these systems 
spin-f and spin- 4- electrons flow along the same physical 
channel; the gate voltage applies in the same way to 
all p aa (a = R/L and a =t,|), so that it effectively 
couples to the charge density p c (x) only. In contrast, in 
TIs the space separation between boundaries allows for 
the presence of two different gate voltages [see Eq.([5)], 
and the helical nature of the edge states leads to the 
coupling to the spin current I s (x) as well. As we shall 
see, this additional handle offers interesting possibilities 
to control the spin current. 



The fact that only p c and I s appear in Eq.(21 ) is due to 
the time-reversal symmetry of the system. Indeed elec- 
tron charge and spin current operators are the only two 
combinations of the four densities p aa that are even un- 
der time reversal, whereas the electron current 

I c (x) = ev F (pnf(x) - PLt(x) + Pri(x) - Pl\{x)) (24) 

and spin density 

Ps{x) = e (pr^(x) + pLt(x) - Pri(x) - Pli(x)) (25) 

are odd. This enables us to highlight the differ- 
ences with respect to the more traditional magnetic 
Aharonov-Bohm effect investigat ed in semiconductors^ 
and graphcnc 7 rings, QHE systems^E±] anc j a i so XI edge 



states^ 1 . In these systems a coupling occurs between the 
vector potential and the charge current I c , and the inter- 
ference is driven by the magnetic flux, breaking time re- 
versal symmetry. In contrast, the Aharonov-Bohm phase 
( 20 1 related to the process in Figj2^b) originates from a 
purely electrostatic effect, preserving time-reversal sym- 
metry: the electron momenta in the two 'arms' of the 



ring are different whenever V 9> t — ^g,B 7^ 0. 

On the other hand, the interferometric process in 
Fig{2jb) also differs from the electrostatic Aharonov- 
Bohm effecpSEZl. While in such effect the electron spin 
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plays no role, here the helical properties of TI edge states 
relate the AB phase 4>ab to spin flipping processes, thus 
modifying the spin current, as we shall describe below. 
A distinction can also be pointed out with respect to the 
Aharonov-Casher effect, observed i n var ious mesoscopic 
semiconductor ring nanostructuref^^. In that case 
a spin precession along the ring arms is caused by 
a uniform Rashba spin-orbit coupling, which can be 
controlled via a gate by varying the asymmetry of the 
quantum well structure. This yields a spin-dependent 
phase difference between the two ring armiPSl. In 
contrast, here spin flip processes occur only locally at 
the QPCs, and 4>ab is actually independent of the spin. 

Finally, another important difference with respect to 
electron interferometers realized with nanotubes or semi- 
conductor rings is concerned with the coupling to the 
biasing electrodes. Such interferometers consist of two- 
terminal setups, where spin-f and spin- 4- electrons are in- 
jected from the same electrode and flow along the same 
physical channel. In contrast, TI edge states are geomet- 
rically separated and the four terminals in Fig{T] allow 
for an independent control of the four injected electron 
species. For these reasons the topology of the interferom- 
eter under investigation here is intrinsically much richer 
than the above interferometers. As we shall show, this 
property enables a tunability of the type (charge or spin) , 
the magnitude and, in some cases, also of the sign of the 
currents. 



and the spin currents fulfill the following relations 

1.(3) - J c (3) 
1,(4) = -J c (4) . 



(27) 



Equation (27) points out that the charge currents mea- 



sured in terminals 3 and 4 are always accompanied by a 
spin current, as expected from the peculiar properties of 
TI edge states. The difference in the relative sign between 
charge and spin currents can be easily understood even 
in the absence of the QPCs. Indeed a positive voltage 
bias applied to terminal 2 injects spin-f electrons into the 
scattering region, whereas a positive voltage bias applied 
to terminal I injects spin- J, electrons, leading to paral- 
lel charge currents and counter-flowing spin currents [see 
Eqs.Q and ([24])]. 

In particular, we are interested in two configurations of 
applied voltage biases V% and V%, namely 



(C) 'Charge'-bias 
Configuration 



(S) 'Spin'-bias 
Configuration 





(28) 



(29) 



IV. RESULTS FOR CURRENTS 

We shall now present the results concerning the cur- 
rents in terms of the applied biases and gate voltages. 
The definitions of charge and spin current arc given by 
Eqs.(24) and (23 1, respectively. Using Eq. Q and ex- 
ploiting the Scattering matrix ( |f 0| ), the currents flowing 
in each terminal are easily evaluated in terms of the Fermi 
distributions fj(E) (j = 1, ...4) of the states injected 
from the four leads 43 . Due to the presence of four dis- 
tribution functions (each including bias voltage and tem- 
perature) and to various parameters of the Fabry-Perot 
interferometer, the behavior of the current exhibits an 
extremely rich scenario. Since we are interested in the 
quantum regime, we shall limit our analysis to the zero 
temperature case. Furthermore, we shall focus on config- 
urations where non vanishing voltage biases are applied 
to the electrodes on the left-hand side of the sample (ter- 
minals I and 2), and the currents are measured in the 
terminals on the right-hand side (3 and 4), which are 
assumed to be grounded, i.e. = Ja = feq, with f eq de- 
noting the Fermi distribution of any lead at equilibrium. 
Under these circumstances, the expression for the charge 
current in the i-th lead (i = 3, 4) reduces to 



IS) 



JdEY, 

J A 1 ' 



\S tJ {E)\\f 3 {E) 



feq(E)) (26) 



The labels 'Charge' and 'Spin' associated with these 
bias configurations originate from the overall degree of 
freedom injected into the scattering region, depicted 
in FigjT] as a dashed box. Indeed in configuration (C) 
the amount of spin-^ and spin- J, electrons injected from 
terminals I and 2 is the same, so that effectively only 
the charge degree of freedom is injected, and no net spin. 
In contrast, in configuration (S) the lead I is negatively 
biased, determining a depletion of spin- J, electrons with 
respect to the equilibrium situation. In this case only 
a spin degree of freedom is supplied to the scattering 
region, with no net amount of injected charge. Similar 
configurations have been discussed in Ref. 42 . In inter- 
mediate situations, i.e. when |Vi| 7^ | V2 1 , both charge 
and spin degree of freedom are involved. The injected 
degree of freedom (charge and/or spin) experiences 
scattering events due to the inter-boundary tunneling 
terms, which determine the current actually measured 
in terminals 3 and 4. For each configuration one can 
define charge and spin conductances in the i-th lead as 
the linear response of I c and I s to the bias V 



dl c (s) (i) 



dV 



i = 3,4 



(30) 



.7=1,2 



In the following we shall describe the behavior of G c and 
G s in the two configurations (C) and (S). 
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A. Charge bias configuration (C) 

Let us start by describing the charge bias configura- 
tion (C) Eq.(28). In this case, using the expressions 
(26l-(p7| for the currents, the definition (301, and the 



properties (17) of the Scattering Matrix, one obtains 



G (3) 
G.(3) 
G c (4) 
G s (4) J 



where denote the transmission coefficients of the 



G (C) = -(T 32 +T 31 ) (31) 



setup. Using Eqs.(18l and (15), one obtains 



U-7? 2 ) 2 



h 1 



Tf 2 + 2T 1 a 2 cos < 



>FP 



(32) 



where T{* 2 is the transmission coefficient of each single 
QPC defined in Eqs.(12), and 4>fp is the Fabry-Perot 



phase. Here we have assumed QPC with equal parame- 
ters, i.e. Tf 2 = T\ 2 . The oscillatory behavior of G^ is 
shown in F igj3| for different values of T" 2 . Recalling the 
expression (19) for 4>fp, one can see that the oscillations 
have a period A(eV) = hvp/L in the sum V 9t B + Vg.T- 
The maxima occur at 4>fp = (2m + 1)tt (to any integer), 
i.e. when the two paths of the FP loop [Figj2^a)] interfere 
destructively, favouring direct transmission from termi- 
nal 2 to 3. These maxima are indeed resonances, for the 
tunneling amplitudes of the two QPC are equal; a slight 
difference in the parameters would make resonances be- 
come maxima very close to 1. This type of behavior of 
G^ for the charge bias configuration (C) is quite similar 
to the gate-induced Fabry-Perot oscillations in the linear 
conductance of a single channe l quantum wire^ or a car- 
bon nanotube transistor! 1 I 2 I 44 I At the level of interference 
processes, the relations between such systems and the 
electron loops originating from spin-preserving tunneling 
processes have been pointed out in Section |III A| Here 
we find that, in configuration (C), the analogy extends 
to the conductance behavior as well. Indeed in config- 
uration (C) the electrodes 1 and 2 are characterized by 
the same Fermi distribution, and so are terminal 3 and 
4. This configuration is thus topologically very similar 
to a two-terminal setup realized with carbon nanotubes 
or semiconductor quantum wires, and the dependence of 
G^ on the Fabry-Perot phase 4>fp can thus be straight- 
forwardly understood in terms of such analogy. 

On the other hand, Eq.(32) has been obtained also 
taking into account spin-flipping tunneling terms. In 
view of the discussion in Section |IIIB[ one would thus 
expect a dependence of G^°' on the Aharonov-Bohm 
phase 4>ab — e (Vg,T — V 9t B)L/HvF as well. Indeed 
such dependence does occur in each of the transmission 
coefficients T 31 and Tv) appe arin g in the expression (31) 
for [see Eqs.(|18[) and ||l5[)]. However, the result 



Eq.(32) shows that this is not the case for the sum 
T32 + T31 and for the conductance G^ c \ which are 
independent of 4>ab- To understand this effect, one can 
consider Figj2] and realize that any AB loop process 
increasing the transmission from terminal 1 to terminal 
3 has a partner AB loop process decreasing the trans- 
mission from 2 to 3, characterized by the opposite AB 
phase dependence, and leading to a perfect cancellation 
of the 4>ab dependence in the sum T 3 i + T 32 . Notice 
that, nevertheless, spin-flipping tunneling amplitude 7/ 
do enter Eq. ( 32 ) through the coefficients T" 2 of each 
single QPC. One can thus conclude that, for the charge 
bias configuration (C), spin-flipping processes affect 
quantitatively the transmission of each single QPC, 
but do not lead to any qualitatively visible interference 
effect between the two QPCs in the conductance G^ C K 
We emphasize that this lack of dependence of 4>ab only 
arises because in configuration (C) terminals (1,2) are 
biased with the equal voltages and terminals (3,4) are 
grounded. It is thus not an intrinsic property of the 
Scattering region, but a specific feature of the charge 
bias configuration (C). As we shall see, the behavior is 
different for the spin bias configuration (S). 



1.00 



0.75 



C 



0.50 



0.25 



0.00 




FIG. 3. (Color on-line) Charge bias configuration (C) 
[see Eq.( |28[ )]. The conductance (31 1 is plotted as a func- 
tion of the Fabry-Perot phase 4>fp = e(V g ,T + V 9t B)L/hvF, 
for different values of single QPC 'horizontal' transmission, 
namely Ti 2 = 0.04 (solid curve), 0.15 (dashed curve) and 
0.48 (dashed-dotted curve), corresponding to spin-preserving 
tunneling amplitudes j p — 0.1,0.2,0.4, respectively, and to 
spin-flipping amplitude 7/ = 0.1. The Fabry-Perot interfer- 
ence pattern allows to switch from the 'off' state (minima) to 
the 'on' state (maxima) by operating with the gate voltages. 
In configuration (C) the conductance is independent of the 
Aharonov-Bohm phase (f>FP = e(V g ,T — V g ,B)L/hvF- 
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B. Spin bias configuration (S) 

Let us now consider the spin bias configuration (S) 
Eq.pg) Making again use of Eqs.Q, §27$, Q, 
and (17 1, one obtains in this case 



G c (3) 1 
G s (3) 
-G c (4) 
G s (4) J 



= G (S) ± e (Ta2 _ r ) (33) 



where are the setup transmission coefficients defined 



in (18). Utilizing the parametrization in terms of the 



single QPC transmission (12), one finds 



4Tf 3 Tf 4 cos0 AS 



i + r 1 a 2 + 2r 1 a 2 cos<^FP 



2 I 



1 - 



0- 



> 



: 



(34) 



1.00 



0.65 



0.30 



-0.05 



-0.40 



-2 



-10 12 3 



FI G. 4 . (Color online) Spin bias configuration (S) [see 
Eq.(29i]. The conductance G (s) (in units of e 2 /h) is plot- 
ted as a function of the Fabry-Perot phase </>fp = e.{Vg,T + 
V g ,B)L/hvF and the Aharonov-Bohm phase 0as = e (Vg,T — 
V g ,B)L/hvF- The tunneling amplitudes for spin preserving 
and spin- flipping processes are -y p = 0.4 and 7/ = 0.2. By 
tuning the gate voltages one can control both magnitude and 
sign of the conductance. Negative conductance values origi- 
nate from spin-flipping inter-boundary coupling. 



Remarkably, a comparison between Eq.([33| and Eq.(32) 
allows to realize that, while the conductance G^ c > ob- 
tained for the charge bias configuration (C) only depends 
on the Fabry-Perot phase 4>fp = e (Vg,T + V 9j B)L/hvF, 
the conductance G^ in the spin bias configura- 
tion (S) depends on both cfipp and the Aharonov-Bohm 
4>ab — e (Vg,T — Vg,B)L/hvp. The latter dependence 
leads to novel features in the TI setup Fig[l]with respect 




-0.5 



FIG. 5. (Color on-line) Slice plots of Figj4] The conductance 
(33 1 in the spin bias configuration (S) is plotted as a function 
of the Aharonov-Bohm phase 4>ab , for different values of the 
Fabry-Perot phase 4>fp, namely 4>fp = (solid curve), 7r/2 
(dashed curve) and n (dash-dotted curve). The dependence of 
G (S) 

on 4>ab results in an Aharonov-Bohm modulation of the 
0FP-dependent Fabry-Perot pattern, and regions of negative 
conductance arise. 



to carbon nanotube based interferometers. The behavior 
of G' s ' is shown in Fig.|4]as a function of 4>fp and 4>ab- 
By varying (frFP, the conductance exhibits a Fabry-Perot 
pattern, qualitatively similar to the case of the charge 
bias configuration (C) shown in Fig. [3j However, by 
varying the Aharonov-Bohm phase 4>ab related to 
gate voltage difference, one obtains a modulation the 
Fabry-Perot peaks, in amplitude and sign. One can 
thus even realize a sign reversal of the conductance 
G^ s \ This effect, absent in the conductance G^ c \ is 
highlighted in Fig[5] where slices of Fig|4] at fixed <j>FP 
values are shown, and regions of positive and negative 
conductance are clearly visible. 

In order to understand the physical meaning of the 
sign reversal and negative conductance, we first analyze 
the problem on the point of view of charge currents. By 
setting the voltage biases in configuration (S) [Eq.dggj)] 
two charge currents are injected, one rightwards along 
the Top boundary and one leftwards along the Bottom 
boundary. The presence of QPCs induces inter-boundary 
coupling and modifies these currents. In particular, 
inter-boundary forward tunneling from one boundary 
opposes to the charge current counter-flowing in the 
other boundary, and may cause the complete blocking 
of the current (i.e. G^ = 0) or even its reversal 
with respect to the applied voltage bias, giving rise 
to nega tive conductance values. This is illustrated by 
Eq.(33), where the conductance G^ is expressed as a 
difference between T32 and T31. Negative conductance 
G^ occurs whenever the 'crossing' transmission T31 
(describing inter-boundary forward tunneling) over- 
comes the 'horizontal' transmission T32 in magnitude^. 
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This effect is in principie present also in a system 
of two usual quantum wires in presence of inter-wire 
tunneling. In TI edge states, however, helicity implies 
that forward tunneling processes necessarily involve a 
spin-flip [see Eq.([3])], thus transferring the effect onto 
the spin currents as well. Notice that in the spin bias 
configuration (S) the counter-flowing charge currents 
correspond to parallel spin currents, since the charge 
current injected along the Top boundary is carried by 
an excess of spin-^ right-movers, and the charge current 
injected along the Bottom boundary by a depletion of 
spin- J, right-movers [see Eq.([23j)]. Thus, in terms of spin 
currents, the tunneling term (131) effectively tends to anti- 
align their directions. This pictorial way also explains 
why negative conductance are not observed in configu- 
ration (C), where the injected spin currents are already 
counter- flowing. In fact, the two AB loop processes 
that mutually cancel in configuration (C) (see Section 
IV A), contribute with the same sign in configuration (S). 



To conclude this section, we wish to emphasize the dif- 
ference between our result of negative conductance and 
the one obtained in Ref.|41j for a corner junction in a 
quantum spin Hall system. In that case the possibility 
that a current can flow out from the lead with the lowest 
applied voltage is an effect due to the strong electronic 
interaction. In contrast, our result shows that negative 
conductance values can be obtained also in the absence 
of electron-electron interaction, due to spin-flipping tun- 
neling terms. Furthermore, we observe that in principle 
negative conductance may occur also in setups with one 
single QPC, like the ones considered in Refs. [30] and 
[32]. However, with one single QPC the value of T 31 may 
not be easily tuned: due to the linear spectrum of TI 
edge states, application of a gate does not lead to electron 
confinement, quite similarly to what occurs in graphenc. 
The QPC parameters are thus essentially determined at 
fabrication level. In contrast, the presence of two QPCs 
gives rise to loop processes, opening the possibility to 
modulate T31 by interferometry instead of confinement. 
Indeed the total 'crossing' transmission T 31 depends on 
both the single QPC 'crossing' transmissions Tg x , and 
on the AB phase 4>ab, allowing to control the sign rever- 
sal. Remarkably, sign reversal only occurs as a function 
of 4>ab and not when varying 4>fp- An intuitive argu- 
ment to explain this property is based on the loop tra- 
jectories depicted in Fig[2] as observed in Section HI C 



to lowest order in the tunneling amplitudes, T31 does not 
depend on Fabry-Perot processes. This property holds 
to all orders, though. 



V. DISCUSSION 

The results presented in the previous section show 
that in both the 'Charge'(C) and 'Spin'(S) bias con- 
figurations a tuning of the conductance is possible 



recall that G^ and G^ represent (up to a sign) the 
linear response of charge and spin currents in terminals 3 
and 4 [see Eqs.(31) and ([33])] . In particular, this implies 
that the interferometer Fig[T] allows for an all-electric 
tuning of the spin current, without invoking magnetic 
fields or polarized ferromagnets. As observed in the 
introduction, this feature represents an advantage in 
view of miniaturization with respect to spin-based de- 
vices involving magnetic materials, due to the difficulty 
in realizing magnetic fields at nanoscale resolution. 
Furthermore, due to intrinsic absence of backscattering 
in TI edge states, the interferometer considered here 
exhibits typical values of the conductance maxima of 
the order of the conductance quantum, a relevant aspect 
for applications as well. In spintronics devices based on 
ferromagnet-normal metal hybrid junctions, for instance, 
the high tunnel energy barriers severely limit transistor 
conductance in the 'on' state, reducing the current 
delivery capability per channel. 

We now wish to comment the specific results obtained 
for each configuration. 



In charge bias configuration (C) Eq.p8| the conduc- 
tance G^ exhibits Fabry-Perot oscillations as a function 
of the FP phase (ppp — e {Vg,T + V g .B)L/HvF similar to 
the case of nanotubes and quantum wires. Notably, in 
this configuration charge and spin conductances are inde- 
pendent of the Aharonov-Bohm phase 4>ab , so that the 
additional degree of freedom V Si t — V 3i b provided by the 
edge channel space separation with respect to the case 
of carbon nanotubes is actually ineffective. In this re- 
spect configuration (C) is topologically very similar to 
the case of other two-terminal charge-based electron in- 
terferometers. Indeed Eqs.(31 1 imply that, if one merges 



the currents flowing into terminals 3 and 4, a pure charge 
signal is obtained, i.e. 

G*°* = G c (3) + G c (4) = 2 G (c) (35) 
G*°* = G S (3) + G S (4) = . (36) 

in accordance with the two-terminal setup analogy. 
Nevertheless, because of the properties of TI edge 
states, the four terminal setup offers the advantage of 
performing spin-resolved measurements, since currents 
flowing in terminals 3 and 4 are carried by a specific 
majority spin component. Indeed configuration (C) 
can for instance be exploited to investigate charge and 
spin tunneling currents through one QPC, as shown in 
Ref. 1521 . where relations similar to (31 1 have been found 



between G c and G s even in the presence of interaction. 

A richer scenario is obtained for the spin bias configu- 
ration (S) [see Eq.(29l], intrinsically different from a two- 



terminal setup. In this case the Fabry-Perot oscillations 
of G^- s ' driven by </>fp are modulated by the Aharonov- 
Bohm oscillations driven by 4>ab — e(Vg,T — V gt B)L/hvF- 
In particular a sign reversal of the conductance is possi- 
ble by varying the gate voltage difference. As observed 



operating with the gate voltages V 9i t and V 3! b- We in Section IV B the sign reversal of the conductance is 
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caused by the spin-flipping terms ^ that effectively be- 
have as an anti-aligning coupling for the injected spin 
currents. Furthermore, we point out that in the bias 
configuration (S) the charge conductances in terminals 3 
and 4 have equal magnitude but opposite signs, whereas 
the spin conductances have equal magnitudes and signs. 
By merging the currents flowing in terminals 3 and 4 

GT = G c (3) + G c (4) = (37) 

G* ot = G s (3) + G s (4) = 2 G (s) (38) 

one obtains a pure spin current, a flow of electronic 
angular momentum that is not accompanied by any net 
charge current.^ The realization of pure spin currents 
represents an extremely attractive problem nowadays, 
for it allows to avoid charge-related spurious decoherence 
effects, and possibly opens new perspectives for coherent 
transport and quantum computing. Indeed a number 
of schemes have been proposed to generate pure spin 
currents.^ On the other hand, the absence of any net 
charge current also makes the detection of pure spin 
currents a difficult task, which requires ad hoc methods 
exploiting e.g. optical technique s] 14 * 26 ! or spin-valve 
effect.^ The setup described here, based on TI edge 
states, is relatively versatile in this respect. Exploiting 
TI edge state properties, it allows to generate spin cur- 
rents by setting the voltage biases in configuration (S), 
to measure them with ordinary methods (ampereme- 
ters) in separate terminals 3 and 4, and to adjust the 
conductances to the desired values by operating with the 
gate voltages. Then, by merging the signals of terminals 
3 and 4, a pure spin current is obtained. The ballistic 
behavior yields values of conductances comparable with 
e 2 /h. These aspects are expected to contribute to make 
the use of spin currents a realistic perspective within a 
short time. 

We now provide some typical values of the relevant 
quantities involved in the setup, which can be realized 
in HgTe quantum wells 32 , where multi-terminal trans- 
port measurements have recently been performed^. The 
period characterizing the oscillatory behavior of the con- 
ductances as a function of 4>fp & n d 4>ab: when expressed 
in terms of gate voltages, equals eAV g T/B — h/vpL. 
Assuming a distance between the two QPCs L ~ 1 fim, 
and a Fermi velocity^ vp — 0.5 ■ 10 6 m/s, one obtains 
a value eAV g ^ T / B ~ lmeV. The possibility for conduc- 
tance sign reversal is related to the order of magnitude 
of T31, the 'crossing' transmissions of the setup. The im- 
portant parameter that ultimately determines T31 is the 
value of the sing lc QPC transmission T 3 a x and T$ v We 
emphasize that relatively weak spin-flipping probabilities 
are sufficient to obtain sign reversal. Indeed, the values 
7 P = 0.4 and 7/ = 0.2 used for Figj4] correspond to sin- 
gle QPC 'crossing' transmissions = = 10% and 
T $2 = T ii = T 32 = T ti = 45%. Assuming two equal 



QPCs, one can see from Eq.(34) that the condition to 



observe a sign reversal can be expressed as the inequality 
(2& - Tf 4 ) 2 - 4T 3 a 1 T 1 a 4 < (39) 



in terms of the single QPC transmission probabili- 
ties. The values of T 31 and Tf 4 can be operatively 
determined by biasing only terminal 1 (Vi — V and 
V-2 = V3 = V4 = 0) and measuring the linear response of 
the currents in terminals 3 and 4, respectively. 



Before concluding, a comment about the effect of 
electron-electron interaction is in order. Similarly 
to other one-dimensional systems like carbon nan- 
otubes or semiconductor quantum wires, in TI edge 
states the Coulomb interaction, which is screened at 
short a nd long distan ces, leads to a Luttinger liquid 
behavioJ^ 1 1 42 ' 52 1 57 1 58 1 . For edge states in HgTe quantum 
wells the Luttinger parameter <?, describing the strength 
of the screened Coulomb interaction, has been estimated 
to range from the weakly interacting li mit g ~ to 
the moderate interaction regime g ~ 

O.PEI] down t0 

the strongly interacting case g ~ 0.5^3, depending on 
both geometrical parameters such as the quantum well 
thickness and the conditions of material growth. The 
analysis carried out here has focussed on the regime of 
non-interacting edge states. This has enabled us to de- 
rive a complete solution of the edge state dynamics for 
arbitrary values of tunneling parameters and gate volt- 
ages, which we expect to qualitatively hold for weak in- 
teraction too. In contrast, in the regime of strong in- 
teraction the Scattering Matrix approach is not applica- 
ble and an exact solution is not available. The use of 
other methods (such as Renormalization G roup An aly- 
sis or perturbative treatments) is mandator} * 29 * 52 * 5 ^ , and 
determining the properties of the setup in the strongly 
interacting case represents a demanding task in general. 
For a single QPC or a corner junction, for instance, it has 
been shown that Luttinger liquid signatur es emerge as a 
power-law behavior in transport propertie d 41 * 42 * 52 * 5 -^. For 
the two QPC interferometric setup in Fig[I] one can spec- 
ulate about interaction effects exploiting some analogies 
with recently studied Fabry-Perot interferometers in the 
strongly correlated regime 4 ^. In particular, a modifica- 
tion of the oscillation periods in 4>fp and 4>ab may be ex- 
pected. Also, in the low bias regime the non-interacting 
limit result may qualitatively be preserved by the finite 
distance between the electrodes, while Luttinger liquid 
signatures should emerge at non- linear transport. How- 
ever, we emphasize that the helical nature of edge states 
in Topological Insulators makes them intrinsically dif- 
ferent from other Luttinger liquid realizations, possibly 
causing the emergence of new features. Furthermore, 
strong correlation makes edge states more sensitive to 
local defects suc h as fl uctuations of ion concentration in 
the doping layers®^ or random bonds at the quantum 
well interfaces^ that may even localize the edge states^. 
For these reasons the determination of transport proper- 
ties of setup Fig[l]in the strongly interacting regime is an 
extremely rich and interesting problem, which deserves a 
separate analysis beyond the purpose of the present pa- 
per. 
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VI. CONCLUSIONS 



In conclusion, we have shown that interferometry of 
TI edge states allows to realize a full electrically con- 
trollable charge and spin transistor, without the use of 
magnetic fields. The charge and spin conductances of 
the proposed setup in Fig|T] can be tuned by operating 
with electric gate voltages applied to the Top and Bottom 
regions between the two QPCs, which modify the phase 
of the two interference phenomena characterizing the sys- 
tem: Fabry-Perot-like loops generated by spin preserving 
tunneling processes and Aharonov-Bohm-like loops gen- 
erated by spin flipping tunneling processes. Two volt- 
age bias configurations have been particularly addressed 
[Eqs.(28)-(29)]. The charge bias configuration (C) leads 
to a conductance behavior that is qualitatively similar 
to the Fabry-Perot oscillations observed in carbon nan- 
otubes electron interferometers, with the additional ad- 



vantage of allowing for a measure of spin-resolved conduc- 
tances as well. In the spin bias configuration (S), where 
counter-flowing injected charge currents correspond to 
parallel spin currents, the Fabry-Perot conductance pat- 
tern is modulated by Aharonov-Bohm oscillations. This 
modulation originates from spin-flipping tunneling pro- 
cesses, tending to anti-aling the spin currents and lead- 
ing to a current blocking or reversal with respect to the 
applied bias, yielding negative conductance values. Fur- 
thermore, this configuration also enables to obtain a pure 
spin current. These peculiar features suggest that the 
proposed TI edge state setup may have a relevant inter- 
est for spintronics applications. 
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